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ABSTRACT 



Height series of Doppler observation at the solar limb 
(covering 3800 - 8700 km distance above the pho- 
tosphere) in H a spectral line obtained by big coro- 
nagraph of Abastu mani Astrophysical Observatory 
l|Khutsishvilil 11986ft show the periodic spatial distri- 
bution of Doppler velocities in spicules. We suggest 
that the periodic spatial distribution is caused by 
propagating kink waves in spicule. The wave length 
is found to be ~ 3500 km. Numerical modelling of 
kink wave propagation from the photosphere to ob- 
served heights gives the wave length of kink waves at 
the photosphere to be ~ 1000 km, which indicates 
to the granular origin of the waves. The period of 
waves is estimated to be in the range of 35-70 s. 
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1. INTRODUCTION 



Most of chromospheric radiation in quiet Sun regions 
comes from spicules, which are jet-like chromospheric 
structures observed at the solar limb mainly in H a 
line. Spicules usually show the group behaviour 
and probably are concen trated between supe rgranule 
cells ( see e.g. reviews of lBeckersI l|1972ft and ISterlind 
( 2000)), thus their formation is connected to intense 
magnetic fields. Typical life time of spicules is 5-15 
minutes. Therefore the propagation of MHD waves 
through the chromosphere with shorter period than 
the spicule lifetime can be traced through spicule ob- 
servations. 

More then 30 years ago lNikolskv and Platoval l|1971ft 
reported that spicules oscillate along the limb with 
characteristic period of about 1 minute and veloci- 
ties about 10-15 km/s. The oscillation, which they 
reported, was just periodic transversal displacement 



of spicule axis at one particular height from the pho- 
tosphere. 

If spicules are formed in thin magnetic flux tubes, 
then the peri odic transverse displac e ment of the axis 
observed by iNikolskv and Platoval l|1971ft probably 
is due to the propagation of kink waves. The kink 
waves can be generated in p hotospheric tubes by buf- 
feting of gra nular motions llRobert sl Il981t iHollwed . 
Il98l iSpruid. IToSll ISasan and KalkofenL 11999ft . As 
kink waves cause the displacement of tube axis, then 
their propagation can be revealed either by direct ob- 
servation of the peri odic displacement of spicule a xis 
along the limb as in INikolskv and Platoval l)1971(l or 
by the Doppler shift of spectral line. The later possi- 
bility arises when the kink oscillation occurs in plane 
of observations. 

The periodic spatial distribution of Doppler veloci- 
tie s in spicules has bee n observed almost 20 years ago 
bv lKhutsishviil l|1986ft (hereinafter paper I) . Unfortu- 
nately, neither observers nor theorists paid attention 
to the phenomenon which can be simply explained 
in terms of kink waves. In this paper, we model the 
observations as propagating kink waves excited due 
to the photospheric granulation. 



2. OBSERVATION 

The big 53 cm coronagraph and universal spec- 
trograph of Abastumani Astrophysical Observatory 
have been used to obtain chromospheric H a and D% 
line spectra at different heights (8 heights) from the 
photosphere (see Paper I) . Instrumental spectral res- 
olution is 0.04 A and dispersion is 1 A/mm in H a . 
Observations were carried out at the solar limb as 
height series beginning at 3800 km height from the 
photosphere and upwards in coronal days when spa- 
tial resolution was 1 arc sec (see paper I). The dis- 
tance between neighboring heights was 1 arc sec, thus 
the spatial distance 3800-8700 km above the photo- 
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Therefore we may consider the velocity distribution 
during one height series as pictured at the same time. 

Thus observations show the evidence of transversal 
wave propagation in spicules. The waves can be ei- 
ther Alfven or kink waves. But as spicules are highly 
structured phenomena (their width are ~ 1 arc sec) , 
the waves probably are kink waves. In next section 
we model the observations as propagating kink waves 
in spicules. 



THEORETICAL MODELLING 



Figure 1. Two height series of Doppler velocity in 
two different spicules in H a spectral line are shown. 
Marked dots indicate to the observed heights (note 
that time difference between consecutive height obser- 
vation is Is). The velocity has clear periodic spatial 
distribution, which can be caused by the propagation 
of kink waves. The wave lengths are around 3500 km. 



sphere was covered. Here we only consider the H a 
observations, however D% line spectra are also inter- 
esting and can be studied in future. The exposure 
time was 0.4 s at lower heights and 0.8 at higher ones 
for H a . The total time duration of one height series 
was 7 s. The consequent height series begins imme- 
diately. The total duration of the observation was 
44 min. Therefore at one height we have a continu- 
ous H a spectra with time intervals of 7 s. Approx- 
imately 300 height series was taken. More details 
about the observations can be found in the paper I. 

We analayzed measured spatial distributions of 
Doppler velocities in H a presented in the Paper 
I. Nearly 20% of spicules show periodic spatial 
distributions of Doppler velocities in height series 
l|Khutsishvili et aLll2004l) . The Doppler velocity spa- 
tial distributions in two height series of two different 
spicules in H a spectral line are shown on Fig.l. The 
periodic spatial distributions are clearly seen. So the 
transversal velocity in that particular spicules is pe- 
riodically redistributed with height, which indicates 
to propagating or standing waves. The wave length 
of propagating waves is 



-3500 km. 



(1) 



The maximal time difference between lowest (3800 
km above the photosphere) and highest (8700 km 
above the photosphere) heights is 7 s, which is the 
duration of each height series. Estimation of the 
distance at which the wave can travel during this 
time gives ~ 350 km (for the Alfven speed ~ 50 
km/s in the chromosphere), which is less than spatial 
resolution of the coronagraph (~ arc sec) , thus less 
than distance between consequent observed heights. 



The mechanism of spicule formation is not well un- 
derstood, therefore we do not concern here how 
spicules have been formed (for ex citation mecha- 
nisms see review of ISterlinel l)200fl|) and references 
therein). As expected kink wave period (with phase 
speed of 50 km/s and wave length of 3500 km the 
period is 70 s) is much shorter than the spicule life 
time, then during this short time interval we may 
consider the spicules as just existing stable struc- 
ture. So we may study the wave propagation into 
this stable structure. Hence despite of very dynamic 
behaviour of spicules, the approximation allows to 
gain main properties of wave propagation. 

Spicules can be modelled as thin magnetic flux tubes 
embedded in field free environment, anchored in the 
photosphere and persisted towards the corona. The 
equilibrium state of stratified atmosphere is defined 
by hydro dynamic vertical pressure balance 



dp e /dz = -p e (z)g, 



(2) 



where p e (z), p e (z) are the hydrodynamic pressure 
and the plasma density through the atmosphere and 
g is the gravitational acceleration. The equilibrium 
state of vertical thin magnetic flux tube also requires 
transversal pressure balance at tube boundaries 



P0(z) + B%(z)/8TT = p e (z) 



(3) 



where po(z) is the hydrodynamic pressure inside the 
tube and B (z) is tube magnetic field. 

Propagation of kink waves in vertical magnetic tube 
embedded in field free environment is governed by 
the e quation (jRobertsl Il981i ISpruitl Il981t iRobertsl 

2004) 

d 2 Z 2 <9 2 £ Po - Pe ,.s 

1^2 =c kir^+9 ; -57, (4 

at az 2 - po + p e oz 

where £ is the transverse displacement, po(z) is the 
plasma density inside the tube, 



Cfc = c A 

is the kink speed and 

CA = 



Pu 
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Figure 2. Schematic picture of propagating kink 
waves in spicules. Due to the propagation of kink 
waves the spectral line will be Doppler shifted peri- 
odically in space if the wave velocity is polarised in 
the plane of observation. The Doppler shift will be 
maximal at velocity antinodes (blue and red shifted at 
corresponding antinodes) and zero at velocity nodes. 



is the Alfven speed. 

Kink waves propagating along the magnetic tube 
lead to the oscillation of the tube axis (see Fig. 2). 
This results in the periodic Doppler shift of observed 
spectral line (in our case H a ) when the kink waves 
oscillate in the plane of observations. At the same 
time the Doppler shift will have periodic behaviour 
in height: it will be maximal at velocity antinodes 
(blue and red shifted at corresponding antinodes) 
and tends to ze ro at velocity nod es (for torsional 
Alfven waves see IZaaarashvilil $2003) ). So kink waves 
give the similar behaviour of Doppler shift as the ob- 
servations (Fig.l). Therefore we argue that the prop- 
agating kink waves are the most plausible reason of 
observed spatial periodicity in H a Doppler velocity. 

Spicule density and magnetic field show almost no 
spatial variation at observed heights (3800-8700 km 
above the photosphere) and the stratification also 
can be neglected, then the equation (4) gives the 
approximate dispersion relation: 
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(5) 



where k z is the vertical wave number and u) is the fre- 
quency of kink waves. From the dispersion relation 
we may estimate the period of kink waves using the 
observed wave length and particular kink speed. The 
typical observed wave lengths are of order ~ 3500 
km (see Fig.l). Then for the kink speed being as 
~ 50 — 100 km/s the expected period of kink waves 
is in the range of 



-35 - 70 s. 



(0) 



The cut-off frequency of kink waves due to stratifi- 
cation at the photosphere is 



ul,. 



2/1(2/9 + 1) 



dh 
dz 



(7) 



Figure 3. The comparison of observations (first plot 
on Fig. 1 ) and the numerical modelling. Vertical axis 
shows the height in km and the horizontal axis shows 
the kink wave velocity in km/s. It is seen that the 
computed kink wave is best fitted to the observation 
if its wave length at the photosphere is ~ 1 000 km. 
In numerical simulation the Alfven speed varies from 
10 km/s at the photosphere up to 100 km/s at the 
height 6000 km. The period of waves is 40 s. 



where h is pressure scale height. For an isothermal 
medium the pressure scale height at the photosphere 
is ~ 125 km, which gives (with (3=1) the cut-off pe- 
riod as ~ 660 s. So the expected period of kink waves 
is well below the cut-off value. Thus the kink waves 
with the period of ~ 35-70 s may easily propagate 
upwards. 

Thus estimations show that observed kink waves 
have higher frequencies than the 5-minute oscilla- 
tion in the photosphere. Therefore p-modes can be 
ruled out to be the source of the waves. This imme- 
diately turns to the granulation. The photospheric 
granulation has been s uggeste d as the source of kink 
waves in thin tubes iSpruitl Il981t iHollwed. Il981i 
lHasan and Kalkofeiill99fll) . Mean granular diameter 
is 800-1000 km, therefore the dynamic granular cell 
probably generates the kink waves with wave length 
similar to the diameter (see Fig. 2). But the observed 
wave length of kink waves (3500 km) in higher at- 
mosphere is few times longer than mean granular 
diameter. The discrepancy probably arises due to 
the increase of kink speed with height: granular cells 
generate the kink waves with the wave length com- 
parable to their diameter, but when the wave prop- 
agate upwards the wave length increases due to the 
increasing phase speed. 

In order to test this assumption we should solve the 
equation (4) along whole solar atmosphere from the 
photosphere to observed heights (~ 4000-9000 km). 
However the height distribution of physical quanti- 
ties are not well known. There are only suggested 
values of density and magnetic field in spicules and 
they are also controversial. Thus the difference be- 
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tween plasma densities inside and outside the spicule 
is not well determined. At the photosphcric level 
the thermal equilibrium inside and outside the tube 
i.e. Tq(z) = T e {z) is a good approximation, which 
gives higher density outside the tube than inside (in 
order to held the equilibrium (2)). But in higher 
atmosphere, say at heights > 2000 km, the plasma 
temperature is definitely lower in spicules than in en- 
vironment. This is because the spicules are seen in 
chromospheric spectral lines formed at much lower 
temperature than the surrounding coronal plasma. 
Also the plasma density seems to be much higher in 
spicules than in surrounding coronal plasma. Due to 
this uncertainty of medium parameters, we take the 
height distribution of Alfven speed and the density 
ratio inside and outside spicule po / p e instead of sepa- 
rate height dependence of density and magnetic field. 
We take the value of 10 km/s for the Alfven speed at 
the photosphere, which increases up to ~ 100 km/s 
at the height 6000 km. Also we take 1/3 for the den- 
sity ratio at the photosphere (being higher outside 
the tube), which increases up to 10 3 at the height 
6000 km. Thus both, Alfven speed and density ratio 
are linear functions of height 

c A = 10 + 1.875^ km/s, (8) 

p e 3 Ah 

where h=125 km is pressure scale height at the pho- 
tosphere and z = corresponds to the solar surface. 

Taking these values of Alfven speed and density ra- 
tio we solve numerically the kink wave equation for 
different wave frequency. And then we fitted the nu- 
merical solution to the observed curves of Dopplcr 
velocity. Result of numerical solution fitted to ob- 
served height series is plotted on the Fig. 3. It is seen 
that the numerical solution is well fitted to the ob- 
servation. In this case the wave period is taken as 

-40 s. (10) 

From Figure 3 it is seen that the wave length of kink 
waves at the photospheric level is 

- 1000 km, (11) 

which is comparable to the granular diameter. It 
clearly indicates to the granular origin of the waves. 



4. CONCLUSIONS 

Height series of Doppler observation at the solar limb 
(3800 - 8700 km distance above the photosphere) in 
H a spectral line show the periodic spatial distribu- 
tion of Dopplcr velocities in spicules. We suggest 
that the periodic spatial distributions are caused by 



propagating kink waves in spicules. The wave length 
at observed heights is — 3500 km. Numerical mod- 
elling of kink wave propagation from the photosphere 
to the observed heights is best fitted to the observa- 
tions when the wave period is ~ 40 s. In this case 
the wave length of kink waves at the photosphere is 
— 1000 km, which probably indicates to the granular 
origin of the waves. 

Our suggested scenario is the following: due to the 
nonstationary character of granular motion the kink 
waves with the wave length comparable to the gran- 
ular diameter are generated in thin magnetic flux 
tubes at the granular boundaries. Then the waves 
propagate upwards along the magnetic flux tube to- 
wards the corona (see Fig. 2) and cause the observed 
height variation of Doppler shift in spicules. 
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